home *** CD-ROM | disk | FTP | other *** search
/ APDL Eductation Resources / APDL Eductation Resources.iso / programs / electronic / rlab / !RLaB / doc / p3 < prev    next >
Encoding:
Text File  |  1994-04-02  |  430 b   |  22 lines
  1. t = (0:10:.05)';
  2. x = exp (-0.5*t) .* (cos (2*pi*3*t) + sin (2*pi*7*t));
  3. X = fft (x);
  4. rfile faxis
  5. freq = faxis (X, .05, 3);
  6. mag = abs (X);
  7. rfile angle
  8. phase = angle (X);
  9. pstart (1,2,"xwin");
  10.  
  11. ptitle ("Magnitude of FFT")
  12. xlabel ("Frequency (Hertz)");
  13. plwid(6);
  14. plot ([freq,mag]);
  15.  
  16. ptitle ("Angle (atan2(imag/real)) of FFT")
  17. xlabel ("Frequency (Hertz)");
  18. ylabel ("Angle (radians)");
  19. plwid(6);
  20. plot ([freq,phase]);
  21. plprint ("p3.ps");
  22.